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ABSTRACT 

Supernova (SN) Refsdal is the first multiply-imaged, highly-magnified, and spatially-resolved SN ever observed. The SN exploded in 
a highly-magnified spiral galaxy at z=1.49 behind the Frontier Fields Cluster MACS 1149, and provides a unique opportunity to study 
the environment of SNe at high z. We exploit the time delay between multiple images to determine the properties of the SN and its 
environment, before, during, and after the SN exploded. We use the integral-field spectrograph MUSE on the VLT to simultaneously 
target all observed and model-predicted positions of SN Refsdal. We find Mg ii emission at all positions of SN Refsdal, accompanied 
by weak Fe ii* emission at two positions. The measured ratios of [O ii] to Mg ii emission of 10-20 indicate a high degree of ionization 
with low metallicity. Because the same high degree of ionization is found in all images, and our spatial resolution is too coarse to 
resolve the region of infiuence of SN Refsdal, we conclude that this high degree of ionization has been produced by previous SNe 
or a young and hot stellar population. We find no variability of the [O ii] line over a period of 57 days. This suggests that there is no 
variation in the [O ii] luminosity of the SN over this period, or that the SN has a small contribution to the integrated [O ii] emission 
over the scale resolved by our observations. 

Key words. Galaxies: high redshift, star formation, ISM, evolution - Gravitational lensing: strong - Techniques: spectroscopic 


1. Introduction 


Supernova (SN) Refsdal ( [Kelly et al.|2Q15| ) is the first observed 
lensed SN resolved into multiple images. Being hosted by a spi¬ 
ral galaxy (Spl 149) that is multiply imaged itself, the discov¬ 
ery presents the possibility of studying the environment of SN 
Refsdal at three different epochs, with lensing affording spa- 
tial scales of -0.7 kpc. Preliminary strong-lensing models (e.g. 
Oguri|[20T5{ [Sharon & Johnson||2015[ [Diego et ^[2015[ ) pre¬ 
dicted that we are now witnessing the second appearance of the 
SN. Therefore, we are able to study the environment of the SN 
before, during, and after its explosion. 

Although SN Refsdal is a very interesting object, determin¬ 
ing the properties of its environment could be even more impor¬ 
tant from the point of view of galaxy evolution. Galaxy growth 
is regulated by feedback and outflows, which are ascribed to en¬ 
ergetic processes such as SNe and active galactic nuclei. The 
necessity of feedback is widely acknowledged, but the detailed 
processes that regulate it are still under debate. Feedback of SNe 
becomes more efficient when the surrounding medium has al 
ready been heated and diluted by young and hot stars (e.g[Stin- 


[son et al.|201 3 [[Hopkins et aH2014t[Artale et al.|2015t[Kim & 


[Ostriker[[2015| [Martizzi et al. [[20151 [Vasiliev et al.[[2015] ), and 
observations of the properties of the interstellar medium (ISM) 
surrounding SNe are crucial to confirm this heating, especially 
at high redshift where feedback is strongest. The resolution and 
luminosity required to study the environment of SNe at high red¬ 
shift can currently only be achieved through gravitational lens¬ 
ing. 

Given its recent discovery, the environment of SN Refsdal 
has not yet been studied in det ail, but its host gal axy was the 
focus of several works. Spl 149 ( [Smith et al.|2009a[ ), located be¬ 
hind the Frontier Fields cluster MACS 1149 (z=0.542), is a spiral 
galaxy at 1.49 with many star-forming clumps in i ts arms, and 
a large magnification ( [Zitrin & Broadhur^[2009[ ). [Yuan et aL] 
( [2Q11|[2QT3] ) found a steep metallicity gradient in Spl 149, mea¬ 
suring a low metallicity (Z <0.3Zo) at the position of SN Refsdal. 
Active, widespread star formation (SF) was shown by rest-frame 
UV observations ([Smith et al.|20Q9a[ ), and through H a detection 


(Livermore et al. 


2012 


2015b at a rate of a few M© yr ^. The size 


and age of the clumps increase inwards, with a young population 
of stars at the position of SN Refsdal ( [Adamo et al.|2013| . 

The environment of SN Refsdal can be studied through the 
presence of certain emission lines. In the optical regime. Mg ii is 
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Fig. 1. Rest-frame spectrum of Spll49 at the position SI an d its centre (SC), w hich have been smoothed over 9 pixels for illustrative purposes. 
For comparison, the red and blue composite GMASS spectra ( |Kurk et al.|2013| ) are also shown, which are normalized to the Spll49 spectrum at 
3000 A. To avoid confusion, we added an offset of (-25, 25) xlO"^^ erg s"^ cm"^ A"^ to the central spectrum and the GMASS spectra, respectively. 
The error estimates are shown by the grey line, while the vertical grey bands are located at wavelengths with significant skyline contamination. 
The short black vertical lines show the wavelength of important Fe and Mg lines. 


frequently used, as it is detectable from the ground over a large 
redshift range. Weiner et al.| ( |2009] ) were the first to report Mg ii 


emission lines in galaxies at z > 1, and these lines have been 


20TTI IMartin et al.||2012HHrb et al. 

120121 IKornei et al.||2013l 

Rigby et al.||2014[|Tang et al.||2014l 

Ben Zhu et al.||2015|). The 


origin of Mg ii emission can either be H ii regions ( see [Kinney 


et al. 

1993 

et al. 

20115 


or resonant scatter ing in outflowing gas ( [Prochaska 


at z = 0.69 that was clearly extended further out than the con¬ 
tinuum emission, in support of the latter scenario. However, |Erb| 
et al.| p012| ) found that the Mg ii emission in their sample was 
more similar to H ii regions, supporting the former. 

In this paper, we study the properties of the region surround¬ 
ing SN Refsdal in the highly-magnified, multiply-imaged, spi¬ 
ral galaxy Spll49 at z=lA9 through high-quality integral field 
spectroscopy, at all of its observed and model-predicted multi¬ 
ple image positions. In Sect. we present the observations, in 
Sect. 1^ we show the results of our observations, we discuss the 
implications in Sect.[^ and we show our conclusions in Sect.j^ 


2. Observations 

We observed MACS 1149 with the Multi Unit Spectroscopic Ex¬ 
plorer (MUSE, [Bacon et al.||2012| ) instrument mounted on the 
Very Large Telescope at Paranal observatory. The total obser¬ 
vation time was 6 hours of the Director’s Discretionary Time 
(DDT), and were executed in service mode (PI: C. Grillo; 
ID 294.A-5032). The telescope was centered at 11:49:35.75, 
^=22:23:52.4, such that three images of Spll49 are within the 
field of view (EOV). The observations were carried out in the 
nights of Eebruary 14th 2015 (1 hr), March 21st 2015 (4 hrs), 
and April 12th 2015 (1 hr). Each hour consisted of two exposures 
of 1440 seconds each, such that the total exposure time adds up 
to 17280 seconds or 4.8 hours. The exposures were taken under 
clear and photometric conditions, and the seeing was <1.1" in 
10 out of the 12 exposures. The last hour of the observations in 
March was executed with a significantly worse seeing of ~2". 

The data reduction was performed using the MUSE Data 
Reduction Software version 1.0. The final datacube consists of 
pixels with a spatial scale of 0.2"x0.2" and a spectral scale of 
1.25 A, although we also constructed a datacube with a higher 
spectral sampling of 0.82 A. We measured the spatial full width 
half maximum (EWHM) of the final datacube from a bright star 
in the EOV to be 0.9", while the spectral resolution is ~2.3 A. 
This seeing corresponds to a spatial scale in the source plane of 
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Cl. 

R.A. 

(J2000) 

Deck 

(J2000) 


^Oii 

(10“^^ erg s 

^ Mg II 

“^cm“^) 

SI 

177.39821 

22.39563 

13.5+^-^ 

39.1 + 1.7 

2.6+0.6 

S2 

177.39771 

22.39579 

12.4+1^ 

29.2+1.0 

2.3+1.0 

S3 

177.39737 

22.39554 


20.4+0.9 

1.6+1.1 

S4 

177.39780 

22.39517 

5 7^^-^ 

-2 0 

16.3+1.2 

1.8+1.3 

SX 

177.40012 

22.39670 

4.8+o| 

15.0+1.9 

1.6+0.8 

SY 

177.40382 

22.40205 

4.0^1 

15.2+1.7 

1.5+1.2 

SC 

177.39702 

22.39599 

' •'-'-0 5 

56.8+1.1 

— 

El 

177.39744 

22.39641 

11-7!o8 

25.8+1.5 

— 

F2 

177.39883 

22.39769 


26.7+1.9 

2.9+1.4 

F3 

177.39975 

22.39656 

9-4!o8 

19.6+1.6 

2.1 + 1.0 

F4 

177.40313 

22.40257 


30.5+1.8 

0.2+0.5 


Table 1. Observed (Sl-4) and model-predicted (future: SX, past: SY) 
positions of SN Refsdal, the centre of Spll49 (SC), and 4 star-forming 
clumps for comparison (FI-4). median values of the model-predicted 

magnification factor (Grillo et al. in prep), emission line flux of [O ii], 
emission line flux of Mg ii X 2796 A corrected for stellar absorption 
(see jGuseva et al.|2013| ). The emission line fluxes are deteremined in an 
0.6"-radius aperture, and are not corrected for magnification. 


~500 pc for a maximum magnification of mt/ = 16 or 4 kpc for 
the minimum magnification yu = 2 in Sp 1149. For details on the 
data reduction, see jKarman et al.| ( [2015j ). 

We extracted individual spectra in 0.6"-radius apertures, 
based on the SN coordinates from Hubble Space Tele¬ 
scope (HST) observations and s trong-lensing -modeling predic¬ 
tions (Sl-4, SX, SY, following |Oguri|[2QT5] ). We updated the 
model-predicted positions of SX and S Y according to our refined 
strong-lensing model (Grillo et al. in prep), and show all posi¬ 
tions and magnification factors p in Table We And no change 
in our conclusions when we vary the positions within the uncer¬ 
tainties of our model, or by adopting positions predicted by other 
authors. The spectroscopic observations span approximately the 
rest-frame wavelength range of 1910 A through 3740 A, which 
includes the spectral lines Fe ii. Mg ii, and [O ii]. 


3. Spectral analysis 

We show the rest-frame spectrum at SI in Fig.[2 where we com¬ 
pare the spectrum to that in the centre of the galaxy (SC) and 
a blue-galaxy composite spectrum from the Galaxy Mass As¬ 


sembly ultra-deep Spectroscopic Survey (GMASS Kurk et al. 
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Fig. 2. The two western images (left), and the eastern image (right) of Spll49. The two figures are divided into four panels; top left: HST 
RGB image with white, yellow, and red circles at the positions of SN Refsdal, the centre of the galaxy, and four star-forming clumps of Spll49, 
respectively, top right: Mg ii emission line flux, bottom left: [O ii] emisson-line flux bottom right: equivalent width (EW) of Mg ii All96 A, where 
a negative EW represe nts emission. The c ircles are at the same positions as the circles in the top left panel. The EW have been corrected for stellar 
absorption (0.5 A, see ( 


Guseva et al. 


20131, and, for displaying purposes, the MUSE maps have been smoothed over a box of 3x3 pixels. 


2Q13| ). The slopes of the spectra at SI and SC are very similar, 
and comparable to the composite of GMASS. This is not sur- 
prisi ng, as different studies h ave shown spread-out SF in Spll49 
(e.g. [Livermore et al.|[2Q12| ). The slight difference in slope be¬ 
tween SI and SC is probably caused by contamination of SI 
from the foreground red cluster galaxy which is responsible for 
the Einstein cross. Looking at the individual lines, the strong 
[O ii] emission doublet is clear, and some other lines show im¬ 
portant differences between S1 and the other spectra. While the 
spectrum at SC and the blue composite GMASS spectrum show 
prominent absorption lines, these lines are significantly weaker 
or absent at SI. This is clearest at the Feii A A 2587, 2600 A, the 
Mg II AA 2796, 2803 A, and the Mg lA 2853 A absorption lines. 

While Mg II absorption is not seen at any of the SN positions, 
Mg II emission is visible at all of them, see Figs.[^[^and Table[^ 
Although the majority of previously reported Mg ii emitters also 
show clear Mg ii absorption, some studies found weak or no ac¬ 
companying absorption ([Weiner et al.|2009}|Erb et al.|2012[[Gir| 


[seva et al.|2013t [Rigby et al.|2014| ). In our case, the absent Mg ii 

absorption, the similarly small rest-frame FWHM of ~120 and 
~140 km s“^ for Mgii and [On] respectively, and the negligible 
velocity offset of Mg ii compared to [O ii], indicate that the Mg ii 
emission originates from H ii regions, rather than from resonant 
scattering. 

In Fig.|^ we show maps of [O ii] and Mg ii emission through¬ 
out Spll49 in all three images. While the [On] emission is 
widespread, the Mg ii emission is local. The clearest Mg ii emis¬ 
sion we see is at position S1, while the center of the galaxy shows 
strong absorption of Mg ii. The measured lower fluxes at posi¬ 
tions S2-4 are probably due to a lower p when integrated around 
S2-4, which differs from the values of point-like images in Ta¬ 
ble [2 Since we are interested in the ISM, we correct the equiv¬ 
alent width and flux of Mgii for stellar absorption (0.5 A for 
young stellar populations, for details, see [Guseva et al.|20~T3| . 

In Fig. we see three additional regions, F2, F3, and F4, 
with significant flux from the Mg ii line, which are located too 
far from the SN to be influenced by it. After extracting spec¬ 
tra of these regions, we find that the spectral properties at F2 
and F3 are very similar to those of SI, see Fig. [^for a com¬ 


parison of the Mg II lines. At F4, we find no significant detec¬ 
tion of Mg II, in an aperture of 0.6", see Table However, we 
find a weak emission line accompanied by more significant ab¬ 
sorption when the size of the aperture is reduced. From the HST 
images we can identify F2 and F3 with two other star forming 
clumps in Spll49. Extracting spectra of other SF clumps does 
not reveal more Mg ii emitting regions, which shows that the ion¬ 
ization conditions vary significantly between clumps in Spl 149. 
We note that although previous work showed that the properties 
of clumps, such as metahicity and stellar age, vary significantly, 
these variations cannot explain the observed large differences of 
the line ratios. 

It is worth noticing that we detect Mg ii emission also at po¬ 
sition SX, which is where the SN has not exploded yet due to 
a longer time delay. This indicates that the high degree of ion¬ 
ization is not caused by the SN, but that it was already present 
before the SN went off. This is strengthened by the detection of 
additional SF clumps with Mg ii emission. We use the ratio of 
[O ii] and Mg ii flux to estimate the ionization parameter U, fol¬ 
lowing ([20^. The emission line flux is determined by 

fitting a Gaussian to the spectrum, while the uncertainties are de¬ 
termined by repeating the calculation on 1000 mock spectra. The 
mock spectra are created by applying a scatter randomly picked 
from a normal distribution with a standard deviation equal to 
the error at every spectral element. We measure a ratio of 10-20 
for Sl-4, SX, and SY, where the main differences are due to 
the relatively low S/N and nearby skylines. The models of [Erb[ 
et al.[ ( [20T^ predict an ionization parameter -2 < log(U) < -1 


for such ratios, whic h is higher than that found for H ii regions 
(-3 < log(U) < -2, [Kinney et al.|199^, and a l ow metahicity 
in agreement with the re sults of [Yuan et al.[ ( [2015] ). A metahicity 
even lower than used by [Erb et al.[p012| )7i.e. Z < 0.2Zo, could 
possibly reconcile the measured ratio with the ionization param¬ 
eter found in H ii regions, although we note that our measured 
ratio is already in agreement with the ratios found by [Guseva 
etakldMn] ), i.e. [O ii]/Mg ii=5-20. 


We find weaker and narrower Eeii* AA 2612,2626 A, at SI 
and S2. This non-resonant emission is accompanied by Ee ii ab¬ 
sorption lines that are weaker than observed in the rest of the 
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Fig. 3. Spectra around the Mgii line at all observed and model- 
predicted SN positions (S1-4,SX,SY), and at three star forming clumps 
(FI-3). The black solid lines show the continuum-subtracted spectra 
which have been shifted to their rest frame using the [O ii] doublet to 
determine the redshift. The vertical dashed lines are at the wavelengths 
of the Mg II doublet, the grey solid lines represent the errors, and grey 
vertical bands mark wavelengths with significant sky emission. The axis 
on the top shows the relative velocity of Mg ii A 2796A with respect to 
[On]. 


Pos. 

Feb. 14th 
(10- 

Mar. 21st 
erg s"^cm 

Apr. 12 th 

■h 

SI 

38.2+3.9 

41.7+2.3 

40.5±5.6 

S2 

28.0+4.7 

28.6+2.2 

32.6+6.5 

S3 

22.4+4.1 

22.4+2.6 

29.2+6.3 

S4 

16.3+6.6 

16.3+2.1 

15.6±5.5 


Table 2. Variability of the [On] emission line at the four observed 
positions of SN Refsdal (SI to S4). The emission line fluxes have been 
measured by fitting a Gaussian to the emission. 


galaxy. Due to the proximity of sky lines and lower magnifica¬ 
tion factors, we are unable to confirm or reject the presence of 
these lines at S3, S4, SX, and SY. 

The [O ii] line offers an opportunity to study a possible in¬ 
crease in brightness of the SN. Since our data were taken in 
three separate months, we can compare the [O ii] fluxes at in¬ 
tervals of 35 and 57 days. We estimate the spatial resolution of 
our observations to be ~ 500 - 700 pc at the positions of the 
SN, using the determined FWHM and the magnification factors 
presented in Table This spatial resolution is less than the res¬ 
olution of HST, and much larger than the region of influence of 
the SN. However, the contrast between the SN region and the 
lensing cluster galaxy is much better at the wavelengths of the 
[O ii] line than in the continuum, and if the SN were the domi¬ 
nant source of [O ii] emission in this 500 pc region, we should 
detect the variability. We study the [O ii] line at the three different 
epochs, and And no significant change in the [O ii] luminosity or 
peak position at the individual positions, see Tablej^ The lack of 
variability means either no variation in the SN [O ii] luminosity 
over this period, or a relative contribution of the SN to the [O ii] 
line that is too small to distinguish with our spectrophotometric 
precision. The small variation is not surprising, as only SN with 
ejecta interacting with circumstellar matter are observed to have 
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clear [Q ii] emission in their spectra (e.g. |Filippenko|19^|Smith 
|et al.|2009bl|Pastorello et al.|2Q15| ). 


4. Discussion 


Most other galaxies where Mgii emission is found in absence 
of absor ption are also lensed by foreground c lusters and galax¬ 
ies (e.g. [Pettini et al.|2()T0[|Rigby et al.|20T4) . Similar to these, 
we And that only some regions of Spll49 show Mgii emission, 
while other regions show clear absorption. This means that if the 
complete galaxy is observed with a single slit, the absorption line 
would dominate the spectrum. This highlights the importance of 
spatially resolving galaxies, which is currently only possible at 
z > 1 for gravitationally-lensed and magnified sources. 


The Mg II emission accompanied by a weak FeiiT 2383A 
absorp tion line, is pre dicted by the mo dels of [Prochaska et al.| 
( |2Q11| ) and observed by |Erb et ar] ( |2012| ). While most Fe absorp¬ 


tion lines are coupled to fine structure transitions, FeiiT 2383A 
is only able to fall back to the ground state. This leads to a larger 
degree of emission line Ailing, similar to Mgii. However, this 
also suggests that the other Fe transitions should be clearly vis¬ 
ible in the spectra, and it remains unclear why recombination in 
the HII regions would affect the resonant Fe ii lines differently 
than the resonant Mg ii lines. 

The absence of the Fe ii lines at all positions of SN Refsdal 
indicates that the medium is optically thin, at least for these tran¬ 
sitions, and that this is not caused by the exploding SN. The low 
opacity at the Fe lines agrees with the occurence of Mg ii emis¬ 
sion from HII regions, as that would also require an optically thin 
gas in front of the emitting regions, possibly leaking continuum- 
ionizing-photons. Assuming TMgn < 1 and Z = 0.2Zo limits 
the column density to A^h < 3 x 10^^ cm"^. This provides an¬ 
other indication that the feedback of hot young stars, possibly in 
combination with previous SNe, has removed most of the cold 
gas from the environment before SN Refsdal exploded, without 
forming the shell required for resonant scattering. The early re¬ 
moval of cold gas is in line with predictions of recent feedback 
models, where young and hot stars dilute the ISM before SN 
explosions, and therefore provides a first step towards the confir¬ 
mation of early-stellar feedback. 


5. Summary and conclusions 

We use VLT/MUSE to study the properties of the ISM around 
the multiply-imaged SN Refsdal. We And clear Mg ii and weak 
Fe II* emission, in combination with almost absent Fe ii absorp¬ 
tion. We determine that the Mgii emission is constrained to 
the clump surrounding the SN, and also And it at the model- 
predicted positions in the other two images. This shows that 
these conditions are not caused by the SN, but are rather the con¬ 
sequence of the young stellar population as a whole. By using 
the ratio of the [O ii] and Mg ii flux, we show that the ionization 
parameter is higher than normally seen in Hii regions. This can 
have important consequences for feedback models, as our results 
suggest to also include stellar feedback. 

We use the [O ii] lines to detect possible variation in the lumi¬ 
nosity of the SN, and thereby constrain the light curve. However, 
we do not And a significant variability, suggesting that the SN is 
either not variable on these timescales, or it is not a major con¬ 
tributor to the integrated [O ii] flux at the scale resolved by our 
observations. 
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